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Abstract. The time variation in the water-vapour bands in oxygen-rich Mira variables has been investigated using 
multi-epoch ISO/SWS spectra of four Mira variables in the 2.5-4.0 ^m region. All four stars show H2O bands in 
absorption around minimum in the visual light curve. At maximum, H2O emission features appear in the ~3.5- 
4.0 /im region, while the features at shorter wavelengths remain in absorption. These H2O bands in the 2.5-4.0 /im 
region originate from the extended atmosphere. 

The analysis has been carried out with a disk shape, slab geometry model. The observed H2O bands are repro- 
duced by two layers; a 'hot' layer with an excitation temperature of 2000 K and a 'cool' layer with an excitation 
temperature of 1000-1400 K. The column densities of the 'hot' layer are 6 X 10^"-3 X 10^^ cm , and exceed 
3 X 10^^ cm~^ when the features are observed in emission. The radii of the 'hot' layer (i?hot) are '^l _R, at visual 
minimum and 2 J?« at maximum, where Rt is a radius of background source of the model, in practical, the radius 
of a 3000 K black body. The 'cool' layer has the column density (A''cooi) of 7 x 10^^-5 x 10^^ cm~^, and is located 
at 2.5-4.0 7?*. A'cooi depends on the object rather than the variability phase. 

The time variation of Rhot/R* from 1 to 2 is attributed to the actual variation in the radius of the H2O layer, 
since the variation in -Rhot far exceeds the variation in the 'continuum' stellar radius. A high H2O density shell 
occurs near the surface of the star around minimum, and moves out with the stellar pulsation. This shell gradually 
fades away after maximum, and a new high H2O density shell is formed in the inner region again at the next 
minimum. Due to large optical depth of H2O, the near-infrared variability is dominated by the H2O layer, and 
the L'-band flux correlates with the area of the II2O shell. The infrared molecular bands trace the structure of 
the extended atmosphere and impose appreciable effects on near-infrared light curve of Mira variables. 

Key words, stars: AGB and post-AGB - stars: atmospheres - stars: variables:general - infrared: stars - stars: 
late- type 



1. Introduction 

Asymptotic Giant Branch (AGB) stars are in the late 
stage of the stellar evolution for low and intermedi- 
ate main-sequence mass stars. Generally, AGB stars are 
pulsating variables as represented by Mira variables. 
According to hydrodynamic model atmospheres, the pul- 
sations lift up matter from the stellar surface and ex- 
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tend the atmosphere (e.g. Bowen 1988). Pulsations create 
shocks, causing a step-like structure in the density distri- 
bution as a function of radius (e.g. Fleischer et al. 1992; 
Hofncr et al. |1998 ). The cooling behind the shock is 



efficient (Woitke et al. 1996) and the temperature de 



creases immediately in the post-shock regions (Fleischer 
et al. |1992| ; Hofner et al. |1998D . 

The extended atmosphere is filled with various kinds of 
molecules. The structure of the extended atmosphere can 



be studied using infrared molecular bands. Hinkle (1978) 



and Hinkle & Barnes (1979) analyzed high resolution spec- 
tra of the oxygen-rich Mira, R Leo, and found two veloc- 
ity components in the near- infrared molecular lines of CO, 
OH, and H2O. They concluded that one component is lo- 
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cated near the boundary region of the photosphere and the 
second component is superposed on the first layer above 
the photosphere. 

However, the molecular bands suffer interference from 
molecules in the terrestrial atmosphere. Recent space- 
borne observations enable more comprehensive studies 
of the molecules in the extended atmosphere. Using the 
Short- Wavelength Spectrometer (SWS; de Graauw et 
al. 1996| ) on board the Infrared Space Observatory (ISO; 
Kessler et al. |1996D , Tsuji et al. ( |1997| ) found CO, H2O, 
CO2 and SiO molecules located above the photosphere. 
Markwick & Millar ( |2000| ) indicated that the 2.8 fim spec- 
tra of a Mira variable consists of two H2 O com ponents 
with 950 K and 250 K. Yamamura et al. ( |l999^ ) identi- 
fied SO2 features in the 7 /im region in three oxygen-rich 
Mira variables with an excitation temperature estimated 
to be 600 K. The band was found to be variable, chang- 
ing from emission to absorption. The time scale of the 
SO2 variation was longer than the period of the visual 
variab le phase. Justtanont et al. ( |1998| ) and Ryde et al. 
( |1999D detec ted CO2 bands in the 12-17 /im region. Cami 
et al. ( ^000 ) found these CO2 molecules fill the region 
between 4-400 stellar radii. Not only molecules but also 
fine-structure atomic lines were found in the upper atmo- 
sphere (Aoki et al. |1998 ). 

II2O is one of the most abundant molecules in the 
oxygen-rich atmosphere and is a large opacity source in the 
near-infrared region. Tsuji ( p^78b| ), using low resolution 
spectra obtained by the Kuiper Airborne Observatory, 
suggested that the 5-8 /J,m region is filled with the II2O 
emission arising from the atmosphere above the pho- 
tosphere (or the hydrostatic atmosphere). However, the 
spectral resolution of those data was too low for further 



study. Yamamura et al. (1999b) found emission features 
from water- vapour bands around 3.5-4.0 /im in o Get, 
which was observed at maximum in the visual light curve. 
In contrast, Z Gas, which was observed near minimum, 
showed absorption features at the same wavelengths. They 
analyzed water-vapour spectra from 2.5 to 4.0 /im with a 
simple 'slab' model. The model consists of two molecu- 
lar layers ('hot' layer and 'cool' layer) with independent 
excitation temperatures, column densities, and radii. The 
'hot' layer with an excitation temperature of 2000 K ex- 
tended to ~2i?* in o Get and stayed at ^li?* in Z Gas, 
where i?* is the radius of the background light source rep- 
resenting the star. They surmised that water layers are 
generally more extended at maximum. 

In this paper we examine the time variation in the wa- 
ter vapour bands in the 2.5-4.0 /im region. Yamamura et 
al. (1999b) analyzed two different stars at different phases. 
To investigate whether the difference in the radii of the wa- 
ter layers is related to the phase difference or not, we an- 
alyzed the spectra of four Mira variables observed several 
times with the ISO/SWS. We found periodical variation 
in H2O bands in ISO/SWS spectra. Some features turn 
from absorption to emission during minimum and maxi- 
mum. This variation is explained by the variation in the 
radius of the H2O layer. These H2O molecules are located 



Table 1. The list of target Mira variables. The properties 
are taken from the General Catalogue of Variable Stars 
(GGVS; Kholopov et al. 198S), and are estimated from the 
light curve within two years obtained by AAVSO (Fig. || 
Mattel |1999|). 



Name 


GOVS 






AAVSO 








Vmax 




Period 


^max 


^min 


Period 




[mag] 


[mag] 


[days] 


[mag] 


[mag] 


[days] 


R Aql 


5.5 


12.0 


284.2 


6.5 


10.8 


274 


R Cas 


4.7 


13.5 


430.46 


6.5 


12.0 


444 


T Cep 


5.2 


11.3 


388.14 


6.1 


10.6 


402 


ZCyg 


7.1 


14.7 


263.69 


8.9 


13.4 


265 


Z Cas 


8.5 


15.4 


495.71 








Get 


2.0 


10.1 


331.96 









in the extended atmosphere. We discuss the variation in 
the structure of the extended atmosphere caused by the 
pulsations. 

2. Observational data 

We extracted SWS observations from the ISO Data 
Archive. In Table ^ we summarize the sampled stars and 
their properties, and in Table |^ we show the observa- 
tion journal. These data are collected from three SWS 
programs: AGBSTARS (P.I. T. dc Jong), TIMVAR (P.I. 
T. Onaka), VARLPV (P.I. J. Hron). The data and re- 
sults of o Get and Z Gas reported in Yamamura et al. 



(19991) are also listed. Four stars (R Aql, R Gas, T Gep, 
and Z Gyg) were observed from 5 to 7 times and the ob- 
served time span covered more than one variability pe- 
riod of each star (Onaka et al. 1999 : Loidl et al. 1999a ). 



These data are ideal to examine the time variation in the 
spectra. The spectra were obtained using the full-grating 
scan mode (AOT 01, scan speed is 1, 2, and 3). The wave- 
length coverage was 2.35-45.2 /im. The spectral resolution 
was A/AA — 300-1000, depending on the wavelength and 
the scan speed. The data were reduced using the SWS 



Interactive Analysis package (e.g. Wieprecht et al. 2001). 
The calibration parameters at October 1999 are used for 
the wavelength, detector responsivity, and absolute fiux 
calibrations. The photometric accuracy of the observed 
spectra is 5-7% for br ight sources in band-1, i.e. below 
4 fim (Decin et al. 200C). Small differences in the ffux level 
between different AOT bands were corrected by scaling 
the fiux at each band with respect to the fiux in band- lb 
(2.6-3.0 fiTo). The correction factor is within a few per- 
cent in band-1. The spectra are re-gridded with a con- 
stant wavelength resolution of A/AA = 300, assuming the 
Gaussian wavelength profile at each wavelength grid. The 
re-gridded spectra are oversampled by a factor of five. 

2.1. An overview of the observed spectra 

The observed spectra are shown in Fig. |^. The positions 
of the major molecular bands are indicated. Due to their 
large optical depth (Fig. ||), H2O features dominate both 
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Table 2. Journal of the ISO/SWS observations. The 
phases are estimated from the optical light curves pro- 
vided by the AAVSO (Mattel |1999| ). at L'-band is cal- 
culated from filter transmission of ESO L'-band (effective 
wavelength is 3.771 /im; van der Bliek et al. |1996D . The 
last column is for the ISO program name. 



Name 


Observed 


Scan 


Phase 


at 


Pro. 




date 




speed 




L' [Jy] 




R Aql-1 


Mar. 


18, 


1996 


1 


-0.08 


778.3 


1 


R Aql-2 


Oct. 


01, 


1996 


2 


0.64 


820.4 


1 


R Aql-3 


Mar. 


03, 


1997 


2 


1.21 


1051.1 


1 


R Aql-4 


May 


03, 


1997 


2 


1.42 


811.9 


1 


R Aql-5 


Sep. 


21, 


1997 


2 


1.99 


741.1 


1 


R Cas-1 


July 


21, 


1996 


2 


0.46 


2156.4 


1 


R Cas 


Dec. 


03, 


1996 


2 


0.78 


2372.9 


2 


R Cas-2 


Dec. 


03, 


1996 


2 


0.78 


2322.1 


1 


R Cas-3 


Dec. 


15, 


1996 


2 


0.80 


2324.3 


1 


R Cas-4 


Jan. 


10, 


1997 


2 


0.86 


2367.5 


1 


R Cas-5 


Feb. 


02, 


1997 


2 


0.92 


2532.4 


1 


R Cas-6 


Aug. 


04, 


1997 


2 


1.35 


2849.4 


1 


T Cep-1 


Aug. 


03, 


1997 


1 


0.38 


2258.7 


3 


T Cep-2 


Oct. 


27, 


1996 


2 


0.58 


1909.8 


3 


T Cep-3 


Jan. 


15, 


1997 


2 


0.79 


1832.4 


3 


T Cep-4 


Apr. 


13, 


1997 


2 


1.02 


2127.7 


3 


T Cep-5 


June 


13, 


1997 


2 


1.17 


2528.6 


3 


T Cep-6 


Sep. 


07, 


1997 


2 


1.39 


2148.1 


3 


T Cep-7 


Dec. 


01, 


1997 


2 


1.61 


1964.0 


3 


Z Cyg-1 


Aug. 


05, 


1996 


1 


0.55 


44.8 


3 


Z Cyg-2 


Oct. 


08, 


1996 


2 


0.79 


34.5 


3 


Z Cyg-3 


Nov. 


24, 


1996 


2 


0.97 


55.4 


3 


Z Cyg-4 


Jan. 


24, 


1997 


2 


1.20 


65.8 


3 


Z Cyg-5 


Mar. 


21, 


1997 


2 


1.42 


55.6 


3 


Z Cyg-6 


May 


15, 


1997 


2 


1.63 


48.2 


3 


Get 


Feb. 


09, 


1997 


3 


0.99 


5773.5 


2 


Z Cas 


Feb. 


26, 


1996 


2 


0.50 


148.9 


2 



Program name 1: VARLPV (P. 
T. de Jong); 3: TIMVAR (P.I. 



R Aql 



I. J. Hron); 2: 
T. Onaka). 



AGBSTARS (P.I. 



RCas 




200 400 600 
JD - 2450000 



200 400 600 
JD - 2450000 



Fig. 1. Th e light curves obtained by AAVSO 
(Mattel 1999). The epochs of the ISO observations are 



denoted as dotted lines. 



the global shape and the small features between 2.5 and 
4.0 /xm. CO2 bands are visible at 2.67 /um. OH absorption 



features are seen around 3.5 /xm in R Aql and T Cep, es- 
pecially near maximum, while these features are not clear 
in R Cas and Z Cyg. 

All of the four stars indicate periodical variations. The 
2.7 /im absorption relative to the maximum flux around 
3.5 /im is generally shallower around maximum than min- 
imum. The flux levels return to a comparable level after 
one period, for example, in T Cep ((/>=0.58, 1.39, and 1.61) 
and in Z Cyg ((/)=0.55, and 1.63). 

Representative ISO/SWS spectra in the 3.8 /xm region 
are shown in Fig. ^ There are three conspicuous features 
which show time variation. The appearance of these fea- 
tures is inverse from minimum to maximum. As we will 
discuss in Sect. 3.2.2, the features denoted with 'x' are in 



absorption and the ones with 'y' are in emission. The defi- 
nition of the absorption or emission at 3.88 /zm is opposite 
to the aspect. 




Wavelength [/J.m] 

Fig. 2. ^Hg^O optical depth (r) calculated from the line 
list of Partridge & Schwenke ( |1997D . Parameters are a line 
width of 5kms^^, a column density of 3 x 10^^ cm~^ mul- 
tiplied by isotopic abundance ratio (0.997), and an ex- 
citation temperature of 2000 K. II2O lines continuously 
occupy the spectrum in the 2.5-5.0 /im region. These 
lines are mainly the vibrational transitions of Awi — 1, 
Awa = 2, and Awg = 1 (Table |). (Author note: This 
figure is modified from the original A&A version, due to 
the restriction of file size of astro-ph. The optical depth is 
smoothed with a resolution of R=2000.) 



3. Model analysis 

3.1. Two-layer 'slab' model 

Water spectra in Mira variables are well represented by 
a model consisting of two water layers (Yamamura et 
al. 1999b). In o Cet, one layer contributes to the emission 
features longer than '-~^3.5/im, and the other layer con- 
tributes to the absorption features shorter than ~3.5/im. 
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R Aql (ISO/SWS) R Cas (ISO/SWS) 



6.0 10 



4.0 102 




3 10 



2 10 



2.6 2.8 3 3.2 3.4 3.6 3 
Wavelength [|im] 

T Cep (ISO/SWS) 



3 3.2 3.4 3.6 
Wavelength [|im] 

Z Cyg (ISO/SWS) 




10^ 



^"2! (|,=i.61 



T Cep 1 (<tf-- 
T Cep 2 (il-- 
T Cep 3 {il-- 
T Cep 4 (<Sf-- 
T Cep 5 (<f-- 
T Cep 6 {(])= 
T Cep 7 ((|)= 



:0.38) 
:0.58) 

fl.79) 

:1.02) 
:1.17) 
:1.39) 
:1.61) 



2.6 2.8 



3.2 3.4 3.6 3.8 




Wavelength [|im] 



3 3.2 3.4 3.6 
Wavelength [|im] 



Fig. 3. ISO/SWS spectra. The optical variability phase ((/>), which is estimated from AAVSO light curve (Mattel 1999 ) 
at each observation, is indicated by the labels. The 2.7 fim absorption with respect to the flux maximum around 3.5 /im 
is generally shallower around maximum than minimum. The flux levels return to a comparable level after one period. 
There are two spectra for 0=0.78 (on the same day) for R Cas, and only R Cas in the program of VARLPV is plotted. 



These two layers are called 'hot' layer and 'cool' layer, 
respectively. 



We use the same model as Yamamura et al. ( 1999bD . 
Plane-parallel, disk shape (or circular slabs), uniform 
molecular layers are superposed, and the radiative trans- 
fer through the layers is calculated. A schematic view of 
this model configuration is indicated in Fig. ||. Each layer 
is described by three parameters, the excitation temper- 
ature (Thot or Tcooi), the column density (A'^hot or A^cooi), 
and the radius of the layer (i?hot or -Rcooi)- The radii of 



the layers are given relative to R.^, , the radius of the back- 
ground continuum source. 

The 'hot' layer has an excitation temperature of about 
2000 K. The water spectrum longer than ~3.5 fim is domi- 
nated by this 'hot' layer (seen in Fig. ^, for example). The 
H2O lines in this wavelength region are mainly caused 
by higher excitation levels of Av2 = 2, Avi = 1, and 
Av3 — 1. The upper 'cool' layer has an excitation temper- 
ature of 1000-1400K and is seen below '^3.5 /im, where 
relatively low excitation levels of Avi — 1 and Avs — 1 
dominate. This 'cool' layer becomes optically thick below 
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3.8 10 



3.4 10 



>^ 3.0 10 



2.6 10 



2.2 10 




TCep((/i=1.39) 



RCas (0=1.35) 



T Cep ((/i=0.79) 



R Cas ((^0.86) 



3.70 3.75 3.80 3.85 3.90 3.95 



Wavelength [|im] 



Table 3. H2O bands of main isotope in near-infrared re- 
gion (Herzberg 1945). H2O has three vibrational levels; 
symmetric stretch (wi), bending of the OH bond (1^2), and 
asymmetric stretch (vs). The lower state of vibrational 
levels is 0,0,0 for all bands. 



Band centre 


Upper state 


[cm-] 




Vl 


V2 V3 


1594.8 


6.270 





1 


3151.4 


3.173 





2 


3657.1 


2.734 


1 





3755.8 


2.663 





1 



discontinuous without any interaction, which are unlikely 
in the real atmosphere. 

For the energy level populations of the molecules, local 
thermodynamic equilibrium (LTE) is assumed. A turbu- 
lent velocity of 5kms~^ is assumed for all the molecules 
(see Yamamura et al. 1999b| ). The line list of H2O is taken 
from Partridge & Schwenkc (1997) and the isotope ratio 
from HITRAN 1996 (Rothman et al. |200lD is applied. The 
line lists of OH and CO2, which are used for data-analysis, 
are also taken from HITRAN. 



Fig. 4. Representative ISO/SWS spectra in the 3.8/im 
region. Three conspicuous features, which show time vari- 
ation, are marked as 'x' and 'y'. These features are in- 
verted between the spectra of T Cep (0=1.39) and R Cas 
((^=1.35) and the ones of T C ep (0 =0.79) and R Cas 
(0=0.86). As we discuss in Sect. [3.2.2 , 'x' is an absorption 
phase and 'y' is an emission phase of the feature. The fea- 
ture at 3.88 fj,m is opposite to its aspect. The flux of T Cep 
is scaled by a factor of 1.5. The dots are ISO/SWS data 
sampling points before re-gridding to constant wavelength 
grid. Although there is a scatter in the data presented, we 
could demonstrate that the phases of the features are cor- 
rect and height of the features generally have uncertainty 
of up to 10%. 



^3.5 /im for a column density of 1 x 10^°-1 x 10^^ cm ^ 
(Table 1^). The column densities in Mira variables usually 



Observer 



exceed these values (Yamamura et al. 1999b; summarized 
in Tabic ||). Thus, the features from the 'hot' layer are 
masked by the 'cool' layer below ~3.5/im. 

The 'disk shape' plane-parallel model (Fig. ||) is a sim- 
plification of the spherical geometry. This simplification 
enables us to treat more than one million H2O lines. The 
validity of this disk simplification will be discussed later 
by comparison with spherical models (Sect. |3.4D . 

We note that 'two layer slab model', or even spherically 
geometric two layer model, does not describe the physics 
in the atmosphere of the star. These models are used to fit 
the observed spectra, and to interpret the observational 
results. Two layers are assumed to be independent and 




Cool Layer 
Hot Layer 



Badkgrciund Source 



Fig. 5. Model configuration. Two disk-shape molecular 
layers are superposed on a background source, and the ra- 
diative transfer is solved though the layers. The top panel 
indicates the side view and the bottom panel indicates the 
geometry seen by the observer. 
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Table 4. The approximate H2O column densities where 
optical depth r becomes unity at each wavelength. The 
turbulent velocity is 5kms~^. 



2cx 


iV-r, » -1 


N-r-, „ -1 
"3.5 (j,m — ^ 


Nt^ , -1 


Nto » -1 






[cm-2] 


[cm-2] 


[cm-2] 


2000 


5.5 X 10^' 


1.4 X lO''' 


3.5 X 10^" 


1.3 X 10^" 


1400 


2.4 X 10^^ 


3.4 X 10^^ 


7.6 X 10^" 


1.2 X 10^" 


1200 


4.9 X 10^^ 


5.4 X 10^^ 


1.0 X 10^^ 


1.3 X 10^° 


1000 


1.1 X 10^^ 


1.0 X 10^^ 


1.3 X 10^^ 


1.5 X 10^" 


800 


3.3 X 10^^ 


2.6 X 10^^ 


1.7 X 10^^ 


2.1 X 10^" 


600 


1.6 X 10^'' 


1.1 X 10^^ 


2.1 X 10^^ 


3.7 X 10^" 


400 


1.6 X 10^^ 


1.2 X 10^* 


2.3 X 10^^ 


1.3 X 10^^ 



Table 5. Results derived by Yamamura et al. (19991:) 



Name 



Thot 

[K] 



?cool 

[K] 



loe 



Nhot 



log A''cool 



-Rhot 



-^cool 



Get 2000f 


1400 21.0 




20.5 2.0 2.3 


Z Gas 2000^ 


1200 22.0 




21.0 1.1 1.7 


^ assumed. 








Table 6. The 


ranges of the parameters for analysis. 


Parameter 


Values 


Grid 


Gomments 


T* [K] 


3000 




Fixed 


Thot [K] 


2000 




Fixed 


Tcooi [K] 


1000-1400 


100 




log iVhotfcm"^: 


20.00-23.00 


0.25 




log iVcooiicm"^ 


20.00-23.00 


0.25 




-Rhot [Rt] 


1.0-3.0 


0.2 


^hot<^cool 


-Rcool [R*] 


1.0-3.0 


0.2 






3.0-4.5 


0.5 




3.2. The model parameters 



We calculated more than 50 000 spectra covering a wide 
range of parameters in order to fit the observed spectra, 
using the analysis. The ranges and step sizes of the 
calculated parameters are summarized in Table ^ There 
are seven parameters. T, N, R of the 'hot' and 'cool' lay- 
ers are given independently. The background continuum is 
assumed to be a blackbody of temperature of 3000 K. 
In the following section the effects of each parameter on 
the spectra are described in detail. We emphasize that the 
radius of the hot layer (i?hot) is an essential parameter for 
the variation between absorption and emission at wave- 
lengths longer than ~3.8/xm. The other parameters are 
less important to the features in this wavelength range. 
In the following figures (Figs. |6- 10) the model spectra are 
presented for a star at a distance of 100 pc and with a 
radius i?* of 3 x 10^-^ cm. 



3.2.1. The spectrum of the central star as background 
source 

Generally, the column densities of H2O layers in Mira vari- 



ables (e.g. Yamamura et al. 1999b) give optically thick 
layers at most wavelengths under consideration. In Fig. ^, 



we show calculated H2O spectra with various tempera- 
tures (r*) of the blackbody as a background spectrum 
of II2O, and we examine the resultant dependence of II2O 
bands on the background source. This figure demonstrates 
that the effects of the background source only appear at 
wavelengths longer than ~3.5/L(m where the optical depth 
of II2O becomes small. In the present analysis, we assume 
that the background spectrum is a blackbody with tem- 
perature r*=3000K, since it is difficult to determine this 
parameter. 



2.4x10 



1.2x10 



6.0x10- 



3.0x10-" - 



1.5x10- 




2.5 



3.0 3.5 
Wavelength [/im] 



4.0 



Fig. 6. The effects of on the spectra are shown. 
The thick lines indicate the synthetic spectra with a 
two-layer 'slab' model, and the dotted lines indicate 
the background continuum level for each T*. From 
top to bottom, r*=3000, 2800, 2600, 2400, 2200 K. 
Other parameters are fixed as Thot=2000K, rcooi=1200K, 
i^hot^l.O-R*, i?cooi=1.4i?,, Miot=l X 10^2 cm-2, and 
-^cooi=l X lO^^cm"^. Due to the large optical depth of 
H2O molecules, the difference in background appears only 
above '-~^3.5 /im. 



3.2.2. The radius of the hot layer 

In Fig. 0, the influence of i?hot is shown. The features 
in the region longer than ~3.5/im vary with i?hot- The 
spectrum shorter than ^3.5 /xm is dominated by the ab- 
sorption by the overlaid 'cool' layer and is insensitive to 
the parameters of the hot layer. When i?hot is «1.6i?* 
with Thot=2000K, the emission from the outer part of the 
layer (r > fills up the absorption originating in the in- 
ner part of the layer (r < i?,) and the spectrum becomes 
almost flat at wavelengths longer than 3.5 /im. Emission 
features appear when i?hot becomes larger than «1.6i?*. 
The effects are most prominently seen in the 3.83 /im fea- 
ture in the top panel of Fig. |7[ Also the slope in the 3.8 /im 
region depends on i?hot- The 3.8 /im region is enlarged 
in the bottom panel of Fig. ^. In this plot, not only the 
3.83 /im feature but also other two features at 3.77 /im and 
3.88 /im are marked. These three marked features are in 
emission, when the features arising from the outer disk 
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of i?hot> r > R^, (region B in Fig. ^) exceed the absorp- 
tion arising from the inner disk of r < i?* (region A in 
Fig. |5|); this is marked as 'y' in the figures. As the fea- 
tures arising from r < are dominant, the features are 
inverted and we indicate the features as 'x'. These three 
features are not formed by a single H2O line but by a com- 
bination of a large number of overlapping H2O lines. The 
definition of the emission and absorption features is oppo- 
site for the 3.88 /im feature, as this feature is formed by 
slightly smaller optical depth compared to the neighboring 
regions. 

3.2.3. The excitation temperature 

The observed spectrum of o Get suggests that Thot must be 
above ~1600K; below ^1600 K, the flux level of the syn- 
thesized spectra shortward of 2.5 /.tm is lower than the flux 
levels of the observed spectra. Because the spectrum of the 
optically thick hot H2O layer behaves like a blackbody 
with temperature of Thot, the flux level around 2.5 fXTn de- 
creases considerably below 1600 K. Thermal equilibrium 
chemistry shows that H2O molecules are abundant below 
~2200K in the atmosphere of red giants and that the H2O 
abund ance d ecreases significantly above this temperature 
(Tsuji 1964 ). The hot layer is expected to have an exci- 
tation temperature in the range of 1600 < Thot < 2200 K. 
With the resolution of 300 it is difficult to determine Thot ■ 
Thus, we fix Thot as 2000 K. 

When the hot layer produces emission features, the 
emergent H2O features always appear similar (Fig. |^). We 
calculated model spectra with Thot of 1800 K and 2000 K 
and confirmed that the other parameters were only little 
affected. i?hot becomes systematically larger by 0.2 i?* at 
most (1 parameter step of i?hot) at Thot =1 800 K compared 
to Thot=2000K. 



3.2.4. The column density 

Fig. |9| shows that iVhot changes the spectra moderately. 
The noticeable features around 3.8 /^m (Sect. I3.2.2D are 
seen when iVhot is in the range of 1 x 10^"'^-3 x lO^^m"^. 
Below or above this range these feature are not detectable 
at the current wavelength resolution. 

The wavelengths shorter than 2.7 /im are dominated by 
the lines in the R-branches, while the lines in P-branches 
start at longer wavelengths. The gap between the two 
branches is seen as an 'emission-like feature' at 2.7 /im. 
The column density of the cool layer (iVcooi) determines 



the feature around 2.7 /im (Fig. 



layei 

i). 



T Cep and R Aql, 



the feature. R Gas and Z Cyg, whose features around 
2.7 /im cannot be clearly seen, have higher column den- 
sities than this range. 



3.3. Other molecules in this wavelength range 
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Fig. 7. Upper: The radius of the 'hot' layer (i?hot) with an 
excitation temperature (Thot) of 2000 K influences the fea- 
tures longer than ^^3.5 /im. The solid lines are by two-layer 
model and the dotted lines are only by the hot layer. The 
effects are most prominently seen in the 3.83 /im feature, 
which is either seen in emission or in absorption depend- 
ing on i?hot- In this case, features longer than '^3.5 /im 
change, and become emission when i?hot is larger than 
«1.6i?,. Other parameters are T,=3000K, Thot=2000K, 
Teooi=1200K, 7Vhot=l X 1022 cm^2^ 7Vcooi=l X 1021 cm-2. 
-Rcooi is given as i?cooi^ = ^hot^ + 1^ [Rl]- Lower: The 
spectra of 3.8 /im region are enlarged. Conspicuous three 
features are marked. Note that a feature at 3.88 /im shows 
'emission like' feature when spectra are in absorption (see 
text). 



3.3.1. OH 

There are sharp OH bands visible in some spectra of R Aql 
and T Gep, especially around visual maximum. The Q- 
branch lines at 2.8 /im do not appear strongly in the ob- 
served spectra (Fig.||). This is because H2O with its large 
optical depth masks the OH Q-branch lines in the 2.8 /im 
region. 

In the current analysis, we include OH molecules only 
in the cool layer, and we assume that the excitation tern- 
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2.4x10 - 




3.0x10-^ I . . . . , . . . . , . . . . 

2.5 3.0 3.5 4.0 

Wavelength [fim] 

Fig. 8. The effects of Thot for the case where the hot 
H2O is seen in emission. From top to bottom T/iot=2000, 
1800, 1600 K. Full lines show the spectra with two-layers, 
and dotted lines show the spectra only with the hot layer 
as a background of the 'cool' layer. Other parameters 
areT,=3000K, rcooi=1200K, i?hot=2.2 i?*, i?cooi=2.4i?,, 
7Vhot=l X 10^2 cm^^ iVcooi=l X 10^1 cm-2. 




2.5 3.0 3.5 4.0 

Wavelength [j-Lm] 



Fig. 10. The influence of A^cooi on the spectrum. The solid 
lines are for two-layer model and the dotted line is for the 
hot layer as the background of the cool layer. From top to 
bottom, 7Vcooi=l X 10^0,3 x 10^0,1 x 102i,3 x lO^^l x 
10^2 3 X 10^2 21^-2 other parameters are r»=:3000K, 
rhot=2000K, Teooi-1200K, i?hot=1.8 i?*, i?cooi=2.0 i?„ 
Ni,ot = l X 10^2 cm-2. 



2.4x10' 




3.0x10-^ - 



1.5x10 I , , , , : , , , , : , , , , I 

2.5 3.0 3.5 4.0 

Wavelength [fim] 

Fig. 9. The effects of A^hot on the spectra when hot H2O 
is seen in emission. The solid lines are for two-layer model 
and the dotted lines are only for the hot layer. From top 
to bottom Nh^'S X 10^2,1 x 10^^,3 x lO^^l x 10^1,3 x 
1020,1 X 102"cm-2. Other parameters are T*=3000K, 
Thot=2000K, rcooi-1200K, i?hot=1.8i?„ i?cooi=2.0i?,, 
Ncoo\=l X 1021 cm-2. 



3.3.2. CO2 

12CI6O2 absorption of t he 10 °l-00°0 and 02°l-00°0 vi- 
bration bands (Herzberg |l945|) is seen at 2.7 fj,m. The CO2 
spectrum of Z Cyg (0 = 0.79) is indicated in Fig. |l^. CO2 
molecules are superposed over the two H2O layers. The 
feature at 2.67 /xm is reproduced well by CO2. The pa- 
rameters of the synthetic CO2 spectrum in Fig. [l| are 
an excitation temperature of 800 K, a column density of 
lxlOi^cm~2^ and a radius of 4.0 i?*. These parameters are 
estimates only. We mask the wavelength region of 2.65- 
2.75 fim for the analysis of the H2O fitting and do not 
discuss the CO2 molecules in this paper. 

3.3.3. SiO 

SiO Av—2 bands at 4.1 /im region are present in the SWS 
spectra. The Au=3 vibration levels around 3 /im region are 
not clear. We mask the wavelengths longer than 3.95 fim 
in the analysis in order to avoid the influence of the 
SiO molecules. 



perature and radius for OH molecules are the same as that 
of cool H2O layer so as to reduce the number of the pa- 
rameters. We confirm that the H2O parameters are not 
changed by more than the typical uncertainties by includ- 
ing OH molecules. 

OH molecules can persist up to 4000 K, which is higher 
temperature than for H2O (Tsuji 1964). The assumption 
of the OH excitation temperature in our model is not ideal. 
The OH lines with higher excitation temperature than the 
hot H2O layer are not seen in the model spectra due to 
large opacity of H2O. We ignore OH with higher excitation 
temperature. 



3.4. Spherical effects 

We adopted disk-shape plane-parallel models to calculate 
a number of spectra with various parameter sets for the 
^2 analysis. We examine the parameters derived from the 
'slab' model. In order to validate this simplification, spec- 
tra of spherical models are compared with 'slab' model 
spectra using the same analysis. In the spherical ge- 
ometry, two H2O layers with uniform temperatures are 
adopted. The density gradient in each layer is assumed to 
vary as . An example is indicated in Table ^. The pa- 
rameters are consistent between these two models. The re- 
sults derived from the simple disk model can provide fairly 
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Fig. 11. CO2 features around 2.7 fim. An example is the 
spectrum in Z Cyg-2 (0=0.79). We put one extra CO2 
layer over the cool H2O layer in the 'slab' model. The 
parameters of synthetic CO2 spectrum are estimated to 
be excitation temperature of 800 K, a column density of 
1 X lO^^cm"^, and a radius of 4.0 i?*. The parameters of 



H2O are hsted in Table |. 



Table 7. Validation of a 'slab' model as a simplification 
of spherical model. A spectrum of spherical model is com- 
pared with 'slab' model using the analysis. Rm and 
i?out are the inner boundary and the outer boundary of 
the spherical model, respectively. Thot is 2000 K for both 
models. 





Tcool 


-Rhot 


floool 


JVhot 


iVcool 




[K] 


[R,] 


[R,] 


[cm-2] 


[cm-2] 






fli„ flout 


flin flout 






Spherical 


1400 


1.0 2.2 


2.2 3.0 


3 X 10^^ 


3 X 10^" 


'Slab' 


1400 


2.0 


2.8 


3 X 10^1 


3 X 10^0 



good estimates of the structure of the outer layers. This is 
because the optical depth of the H2O layer is so high (op- 
tical depth at 3.5 /im becomes about unity with a column 
density of 1 x 10^"'^ cm~^ and an excitation temperature of 
2000 K; Table that only the surface of each sphere is 
important for the synthetic spectra. The radius where the 
optical depth becomes thick in a spherical model is repre- 
sented as the radius of the 'slab' model. If we compare the 
spectra with these two geometric models in an optically 
thin case (say column density of 10^^-10"^^ cm~^), we find 
larger discrepancies. 



3.5. Fitting process 

The best fit for each observed spectrum is derived by the 
analysis. The spectral resolution is 300-500 in these 
wavelengths. All the spectra are convolved by a Gaussian 
with the wavelength resolution of A/AA — 300. Xm is the 
reduced x^, and divided by the degrees of freedom m. 
The best parameter sets are found by minimizing Xm- The 
wavelength range 2.5-3.95 /im is used except for the 2.65- 
2.75 ^m region which is strongly affected by CO2 bands. 



The SWS scans the same wavelength region twice. The 
difference between two scans is a good indication of the 
observational uncertainty. Two scans are normally con- 
sistent within 0.6 % of the flux over 2.5-4.0 iim. In some 
observations AOT bands ID and IE (3.5-4.0 /zm) show a 
significant difference between two scans. The difference is 
as large as 10 % of the flux in the worst case (R Aql-1; we 
discard this observation from the later discussions). We 
added this difference in flux and 1 % of the flux, which 
represents the scatter of the data around the re-gridded 
data, to an observational error, cr. The latter part is neces- 
sary to avoid diverging to the infinity. The error of each 
parameter is obtained by calculating grid points within 3 
a level from the parameter of the best x^, when other pa- 
rameters are fixed. The error of our fit is not a continuous 
function but rather a grid scale. 



4. Results 

4.1. Tlie best-fit spectra and their parameters 

The best fit synthetic spectra corresponding to the ob- 
served spectra (Fig. ^ are shown in Fig. [ij. Our simple 
'slab' model reproduces the global shape of the observed 
spectra. Residuals between the observed spectra and the 
fitted spectra are shown in Fig. |l^. The fit is as good as 
a few % in the most of the investigated spectral range 
and about 10% at the worst part. The 3.8 ^m region of 
the observed spectra and the model spectra is enlarged in 
Fig. ^ The variation with phase in the three conspicuous 
features is seen in the observed spectra, and these features 
are reproduced by the model except at (j)—0.78 and 0.80 
of R Gas (R Gas, R Gas-2, R Gas-3). The small features of 
the observed spectra are generally reproduced in respect 
to the trend of emission or absorption. However, the de- 
tailed shapes of the small features (e.g. height or width) 
are sometimes not fitted well. This implies that two layer 
slab model is too simple approximation and the atmo- 
spheric structure is more complicated in Mira variables. 

The parameters of the best fit spectra are listed in 
Table ||. All four stars show the common tendency that the 
radius of the hot layer (i?hot) is about 1.0 i?» around min- 
imum and about 2.0 i?, around maximum ( Fig. |l4|) , con- 
firming the suggestion in Yamamura et al. ( 1999b|) . The 
cool layer is located at about 2.5-4 i?*. The column den- 
sity of the cool layer (A^cooi) varies from object to object 
rather than showing a trend with the phase and log A^cooi 
ranges from 20.75 cm~^ to 22.25 cm~^. 

The column density of the hot layer (A'hot) has a large 
uncertainty as discussed in Sect. [3.2.2 . From the equiv- 
alent width, we could give a lower limit when the fea- 
tures are clearly seen as emission. At the phases with the 
strongest emission features at 3.83 /im in each star (R Aql- 
5 at <?!)=1.99; R Gas-5 at (/)=0.92; T Gep-3 at (/)=0.79; 
Z Gyg-3 at 0=0.97), the column density of the hot layer 
should be higher than «10^^ cm~^, otherwise the features 
would not have been detected. 
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Fig. 12. The fitted results and the residuals. Two spectra at (f)—0.78 (on the same day) and one spectra at 
R Cas are fitted with the same parameters, and the synthetic spectra are identical. 



=0.80 of 



There is an uncertainty remaining in the fits of indi- 
vidual spectra. The spectra of R Cas, R Cas-2, R Cas-3, 
which were taken within 2 weeks, are fitted by the same 
parameter sets. These three synthetic spectra based on 

analysis show absorption features in the 3.8 /zm region. 
However, all of them show emission features in the ob- 
served spectra (Fig. |l^). Thus, the fitted parameters for 
these three spectra are not reliable. Also the features in 
the 3.8 fim region in Z Cyg-6 are not reproduced with the 



best parameters. However, the data of Z Cyg-6 are not af- 
fected by the instrumental error; the features in the 3.8 /zm 
region in Z Cyg-6 seem to be real. Thus, the parameters 
of Z Cyg-6 are doubtful. 

In some phases, the features in the 3.8 /im region do 
not clearly appear in the observed spectra, for example, 
in T Cep-2 (0=0.58). There are at least three solutions: 
(1) The column density of the hot layer is extremely high, 
and the H2O lines are completely saturated. (2) The radius 
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layer is so low that the features are undetectable. We can- 
not distinguish the three cases from the observed spectra. 
Although the analysis suggests that case (1) is the most 
likely for most of such spectra, cases (2) and (3) cannot 
not be ruled out. 

4.2. The 3.8 fim region 

The most interesting parameter in this study is the radius 
of the hot layer, i?hot- The observed spectra in Fig. |l^ 
indicate that the 3.8 /im features show emission around 
maximum and absorption around minimum. This suggests 
that i?hot is larger than wl.Gi?* around maximum and 
smaller than that around minimum. The features in emis- 
sion phase are not always exactly the inverse shape of the 
absorption features, possibly because of the complicated 
structure of the atmosphere in Mira variables. 

The variation in i?hot derived from the model fitting 
is plotted in Fig. |l^ against the optical phase. In all four 
cases, i?hot reaches about 2 i?* around maximum and be- 
comes ^1 i?, around minimum. 

There are slight differences among the stars. Z Cyg still 
shows absorption features at phase 0.79, while the other 
three stars already show weak emission. This trend ap- 
pears to be related to the characteristics of the light curve 
(Fig. I). The hght curve of Z Cyg shows a steep increase 
and a gradual decrease in the V-band. The visual magni- 
tude of Z Cyg at 0=0.79 is near the minimum, while in 
other stars the visual magnitudes already start increasing 
at this phase. Fig. ^ shows -Rhot as a function of visual 
magnitude at each epoch, after subtraction of the maxi- 
mum magnitude. A separation is found at V — Vmax of 
3.0 mag; below that i?hot is ~2i?* and above that i?hot 
is 1.0 i?*, indicating the variation in i?hot is connected 
with variability. One of the classification criteria of Mira 
variables is that visual amplitude is larger than 2.5 mag. 
Fig. ^ implies that the emission features might be related 
with the large-amplitude pulsating variables. 



Fig. 13. The spectra in the 3.8 /xm region are enlarged; 
the ISO/SWS spectra (bold hues), and the model spectra 



(thin line). The three conspicuous features (see Sect. 3.2.2; 
Fig. which shows time variation, are marked. The 
mark, 'x' shows the absorption phase and 'y' shows the 
emission phase of ISO/SWS spectra. The definition of 
the emission and absorption feature at 3.88 /xm is oppo- 
site to their aspects. Around minimum the features are 
observed in absorption, while around maximum they are 
seen in emission. The error bars below the phase indicate 
a (Sect. |3^ ) at 3.83 /xm. Shght wavelength shift of the fea- 
ture peak is seen, especially in Z Cyg. These shifts seem 
to be observational errors. 

of the hot layer just falls at «1.6i?* where the contribu- 
tions of emission fills up the absorption features. (3) The 
hot layer disappears, i.e. the column density of the hot 



5. Discussion 

5.1. The dynamical motion in the extended 
atmosphere 

The spectra of II2O bands in the 3.8 /im region are re- 
produced well assuming rhot=2000K. That temperature 
indicates that the hot layer traces a relatively inner part of 
the extended atmosphere. The variation in i?hot indicates 
the variability of the structure in this region. 

The variation in i?hot can be interpreted by the varia- 
tion in the density structure in the extended atmosphere 
due to pulsations. Pulsations create local high-density 
shells or shocks in the extended atmosphere. If the gas 
temperature in the shock is high enough to dissociate II2O 
molecules, they are re-formed behind the shock front: in 
this case i?hot traces the region with the highest H2 O den- 



sity located behind the shock front. Duari et al. (1999) 



indicated that II2O can be reformed in a post-shock re- 
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Table 8. The best fit parameters, and the estimated uncertainties. OH is included in some fits with a fixed column 
density of 1 x 10^^ cm^^. The numbers of errors are not given when the errors are smaller than 1 step of the model 
parameters. 



Name 


Phase 


^cool 


-Rhot 


-Rcool 


log 


' iVhot 


log 




log iVoH 


Xm 






[K] 


[R.] 




[cm-2] 


[cm-2] 


[cm-2] 




R Aql-1 


-0.08 


1400 




7^ — .1(1 ') — 

^•^-0.2 


21 


7C+U.5U 
— 0.50 


21 


oc + i.UU 
— 0.50 




1.1 


R Aql-2 


0.64 


1400 


2.2 


2 5+0. 2 


22 


50+°-^? 


21 


05 + 0.25 
^—0.25 




2.5 


R Aql-3 


1.21 


1400 


2 Q+0.2 


3.0-0.2 


22 


9c;+0.75 


20 


-75+0.50 
—0.50 


21 


2.2 


R Aql-4 


1.42 


1200 


1 0-^0.2 
'-■'^-0.2 


2.4 


21 


50+°-?5 


21 


'^'-'—0.50 




2.4 


R Aql-5 


1.99 


1400 


^•"-0.2 


2 5+0. 2 


21 


75+0-^^ 

— U . z o 


21 


— U.Zo 


21 


1.5 


R Cas-1 


0.46 


1100 


1.2 


2.4 


21 


' — 0.50 


22 


— \j.2o 




2.7 


R Cas 


0.78 


1000 


1.0+ 


3.5 


20 


yr+0.25 

— 0.75 


22 


— U.ZO 




2.9 


R Cas-2 


0.78 


1000 


1.0+'''' 


3.5 


20 


-TC+O.SO 
— 0.75 


22 


— \J.2o 




3.1 


R Cas-3 


0.80 


1000 


1 n4-0 6 

1.0+ 


3.5 


20 


'-^ — 0.75 


22 


00+^f, 




3.2 


R Cas-4 


0.86 


1200 


2.2 


3.5 


22 


25+0.50 
— 0.25 


21 


7^+0.50 
^—0.25 




2.5 


R Cas-5 


0.92 


1100 


2 2+0-2 


4.5 


22 


25+0.75 

— 0.25 


21 


50+°il 

— U.Zo 




2.3 


R Cas-6 


1.35 


1100 


1.2^ 


2.6 


21 


— 0.75 


21 


— \J.2o 




3.1 


T Cep-1 


0.38 


1200 


1.0+ 


2.4 


21 


— 0.50 


21 


50+r.n 

— 0.50 




4.1 


T Cep-2 


0.58 


1300 


1.8 


2.6 


22 


— O.oO 


21 


50+"1° 

— O.zo 




2.2 


T Cep-3 


0.79 


1300 




2.6 


21 


^^ — 0.50 


21 


"" — 0.^5 


21 


2.3 


T Cep-4 


1.02 


1300 


1.8 


3.0 


22 


— V.zo 


21 




21 


2.5 


T Cep-5 


1.17 


1400 


1.8 


2.4 


22 


— 0.25 


21 


25+0.50 

—0.25 


21 


3.4 


T Cep-6 


1.39 


1200 


i-2io:2 


2.2 


21 


oc:+0.25 
^^-1.00 


21 


c;n+0.50 
>^"-0.25 




3.7 


T Cep-7 


1.61 


1000 


1.0+" 


3.5 


20 
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gion. If the shock is moderate and H2O molecules do not 
dissociate, the H2O density increases at the shock front 
and i?hot traces the motion of the shock front. In either 
case, the variation in i?hot traces the motion of shocks. 

Around minimum, a high H2O density shell is created 
by the pulsation shock near the stellar surface; the area 
of the extended region is smaller, and the contribution of 
the emission component originating from the extended re- 
gion is weaker than the absorption component (Fig. |l6|) . 
Toward maximum, the high II2O density shell moves away 
from the stellar surface. The hot layer is extended and 
emission is observed. This high density shell continues to 
expand after the maximum, but its temperature and den- 
sity decrease. The contribution of this shell to the features 
in the 3.8 /im region gradually decreases. At the next min- 
imum, a new shock is formed on the surface and produces 
the absorption features. 

The radius of the hot layer is measured in units of i?, ; 
the radius of the background continuum source (in prac- 
tice, a 3000 K blackbody). It is uncertain how the stel- 
lar radius varies with phase. Interferometric observations 



show that the stellar radius of Mira variables is rather 
stable, and varies only by less than 7% at 902 nm, a 
wavelength little contaminated with absorption (Tuthill 
et al. 1997 ). The Rossela nd ra dius derived from theoreti- 
cal models (Bessell et al. 1996 ) varies by about 20%; the 
radius is smallest at maximum. These variations in the 
stellar radius are smaller than our derived variations in 
Rhot/R*, which varies by a factor of two. Thus, the loca- 
tion of the hot H2O layer actually varies by as much as a 
factor of 1.5-2.0. 



Hydrodynamic model calculations indicated that II2O 



molecules are in non-LTE (Woitke et al. 1999). We can- 



not discard non-LTE excitation for the cause of emission 



features. Recently, Schweitzer et al. (200C) developed a 
method to solve non-LTE for H2O molecules. Such an ef- 
fect should be considered for future work. 
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Fig. 14. The time variation in the radius of the hot layer 
(^hot)- When _Rhot is larger than 1.6 i?,, the spectrum 
around 3.8 /im is in emission, while i?hot smaller than 
1.6 i?* indicate absorption features. All four stars indicate 
about 2.0 i?* around maximum and about 1.0 i?* around 
minimum. 
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Fig. 15. _Rhot as a function of visual magnitude at each 
epoch subtracted from maximum magnitude. There is a 
separation; when V — Vmax is larger than 3.0 mag -Rhot is 
'^l.Oi?*, and smaller than 3.0mag it is ^^2.0 7?*. 



5.1.1. Comparison with the structural variation in the 
extended atmosphere traced by other molecules 



Hinkle et al. (|1982| [1984D and Lebzelter et al. ( [1999D found 
periodical variations in the near-infrared CO lines in the 
H-band and K-band of several Mira variables. Hinkle et 



al. ( |1982| , |198J) found that the CO Aw = 3 lines, with ex- 
citation temperatures of 2200-4000 K, show a P Cyg pro- 



Hot 



Around minimum 
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Temperature and 
density are decreasing 



1 R:i 



2 R* Radial direction 



Fig. 16. A schematic view of the motion of the high 
H2O density shell. From top to bottom, it is in the time- 
sequence. Around minimum, a high-density shell is located 
on the surface of the star and H2O is observed in absorp- 
tion. The high density shell is moving outwards, and lo- 
cated at «2 i?* around maximum. The emission features 
arise from the high density region expanded further from 
the 'photosphere'. After the maximum, the high density 
shell continues to go outwards. As the shell is expanded, 
the temperature and density decreases and its contribu- 
tion to the spectrum in 3.8 /im region is getting lower. 

file; the blue-shifted emission feature is visible from (p^O.G 
to 0.8. The CO column density is lowest around 0=0.85. 
These variations are clearly different from those of hot 
H2O. The dissociation energies of CO and H2 O mo lecules 
are 11.09 and 5.11 eV, respectively (Dappen 199£ ). It is 
reasonable to assume that the CO molecules in Hinkle et 



al. (1982, 1984) trace even a hotter region interior to the 
hot H2O layers. The authors explained the behavior of 
the CO lines with the shocks caused by the pulsation. A 
shock affects the CO layer around (/)=0.7. The high den- 
sity region caused by the shock influences the H2O layer 
around maximum. It is qualitatively reasonable that this 
phase lag of the expected shock penetration appears be- 
tween the CO and H2O layers because of the difference in 
transition energy and binding energy of these molecules. 

Yamamura et al. (1999a) found SO2 in oxygen-rich 
Mira variables and examined the time variation. Unlike 
H2O, the time variation in the SO2 band does not corre- 
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late with the visual variability but changes from emission 
to absorption over more than a single period. The location 
of the SO2 layer is estimated to be about 4-5 i?, in T Cep. 
It is implied that the inner region in the extended atmo- 
sphere (hot H2O layer) follows the optical variability, but 
the outer region (SO2 layer) does not. It cannot be ruled 
out that the outer region has a longer variability period 
than the inner region because the pulsation shocks do not 
directly influence the outer region. 



5.1.2. Comparison with the hydrodynamic models 



higher than w 1 x 10^^ cm^^ to detect the emission features 
in the 3.8 /im region. The density calculated by Bessell et 
al. (1996) exceeds the column density 1 x 10^^ cm~^ out 
to 0.1 stellar radii. Here we assume 3 x 10^^ cm for one 



Bessell et al. ( |1996D calculated atmospheres of oxygen-rich 
Mira variables including hydrodynamics. They show CO 
line profiles for several transitions. The CO first-overtone 
(excitation level of 0.2 eV) is seen in emission at 0=0.8- 
1.0. They explained that the CO emission was produced 
in the pre- and post-shock region in the in ner part of the 
extended atmosphere. Winters et al. ( 2000 ) calculated the 
CO line profile for both Av — 1 and Av = 2 transitions. 
CO lines, especially those of Aw = 1, show a P-Cyg profile 
with a strong emission component. Although their model 
is calculated for carbon stars, the general trend may be 
applied to oxygen-rich stars. These theoretical models in- 
dicate that CO lines are seen in emission due to the high 
density regions of CO molecules lifted by the pulsation 
shocks. The appearance of H2O lines can be interpreted 
in a similar manner. 



Bessell et al. (1996) showed the temperature and den- 
sity structures in the extended atmosphere. A shock oc- 
curs at the photosphere and moves outwards. Only mod- 
els for the fundamental mode of pulsation can produce 
shocks in the extended atmosphere. Their models indi- 
cated that the temperature does not increase drastically 
at the shock front (nearly isothermal shock). According to 
their Z-modcl (fundamental pulsation mode and luminos- 
ity of 6310 Lq), the location of the shock increased by a 
factor of 1.8 from minimum to maximum. The radius of 
the II2O layer measured in our work is a factor of 1.5-2.0 
larger around maximum than around minimum. Assuming 
that the hot II2O layer traces either the location of shock 
itself or post-shock region, the approximate scales of the 
shock propagation are comparable with that model. 



Bessell et al. ( 1996 ) also calculated the temperature 
and the density around the shock region. The temper- 
ature and the density at the shocked region have uncer- 
tainties according to Bessell et al. ( 1996 ). Here, we discuss 
the post-shock region. They show one example, D-model 
(fundamental pulsation mode and luminosity of 3470 Lq). 
According to one of their D-models, the gas temperature 
is 2000 K and the mass density is 1 x 10~^^ gcm"'^ at the 
post-shock region at maximum phase. The gas tempera- 
ture is comparable to the excitation temperature of our 
hot layer. The H2O number density becomes 5 x 10® cm~'^ 
if we assume the gas is dominated by H2 and the abun- 
dance ratio of H2 and H2O is 10~^. Our analysis indicates 
that the column density of the hot layer is required to be 



stellar radius. Thus, the density in Bessell et al. (1996) is 
consistent with our estimate. 

Another hydrodynamic model is presented by Woitke 
et al. (1999). Unfortunately, the phase is not given in 



their work, but a shock with a gas temperature of 2000 K 
reaches 2 stellar radii. In their model, the H2O num- 
ber density at the shock front is only 3 x 10^ cm^"^. A 
shell thickness larger than 10^^ cm (more than ~30 stellar 
radii) is required to reach the H2O column dens ity o f 
1 X 10^^ cm^^. The number density in Woitke et al. (1999) 
is by more than an order of magnitude lower than our re- 
sults. Their results may be applicable to the observations 
of irregular variables and semi-regular variables by Tsuji 
et al. ( 1997 ), and not to stronger pulsating stars such as 
Mira variables. 



5.2. The optically thick H2O layers in oxygen-rich Mira 
variables and their infrared variability 

The near-infrared variability has been measured for Mira 
variables. For example, the L'-band magnitude of R Aqr 
varies about 0.6 mag (Le Bertre 1993| ). In our sample, the 
variations are about a factor of 1.4-2.2 (0.5-0.7 in magni- 
tudes) at L'-band (Table ||). As the opacity is dominated 
by H2O in this wavelength range, it is likely that H2O 
molecules play an essential role in the flux variation. 

From Figs. |€ - 10 , the parameters which drastically vary 
the absolute flux at 3.8 fim (t he effe ctive wavelength of the 
L'-band; van der BHek et al. |1996D , are the radius of the 
hot layer (i?hot; see Fig. 0) and the excitation tempera- 
ture of the hot layer (Thot; see Fig. ||). The variation in 
i?hot from 1.0 to 2.0 i?* induces a variation in the emis- 
sion by a factor of 4.0. The difference in the flux between 
Thot = 1600 and 2000 K is a factor of two. The flux varia- 
tion is explained by the larger radius and lower excitation 
temperature of the high H2O density shell at maximum 
compared to minimum. 

Tuthill et al. (2000) measured the radius of a Mira 



variable with interferometry: they mentioned that the ra- 
dius at 3.1 fim is larger than those in other near-infrared 
bands, due to molecular opacity. The hot H2O layer in our 
work fits this opacity source. 

These optically thick H2O layers are important for ra- 
diative transfer in the extended atmosphere. Model cal- 
culations of the extended atmosphere with precise treat- 
ment of molecular opacity (e.g. Hofner 1999 ; Helling 2000) 
must be compared with the observations to obtain better 
understanding of the physical conditions of the extended 
atmosphere in pulsating stars. 

6. Conclusions 

We reported the time variation in H2O bands in four 
oxygen-rich Mira variables. The near infrared water- 
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vapour bands at 2.5-3.95 /Ltm follow the periodical vari- 
ation. Emission features are seen at ~3.5-3.95/uni around 
maximum while absorption features are detected around 
minimum. The spectra are well fitted with 'slab' models, 
which consist of two H2O layers (hot layer and cool layer). 
The radius of the hot layer varies from ^li?* to ~2i?* 
during visual minimum and maximum. 

The periodical variation in the features arising from 
the extended atmosphere suggests that the structure of 
the outer atmosphere is varying with the pulsation. The 
pulsation produces a shock in the atmosphere, and the 
hot H2O layer traces the high H2O density caused by the 
shock. The high H2O density shell expands from inside of 
the extended atmosphere to outwards from minimum to 
maximum. In our analysis, i?* is measured as the relative 
number of the stellar radius. Considering the variation in 
the radius of the star, the high density shell is by a factor 
of 1.5-2.0 more extended at maximum than at minimum. 

Due to large optical depth of H2O bands, the spectra 
in the 2.5-4.0 /itm region are dominated by the H2O in the 
extended atmosphere. The flux variation in the L'-band 
is primarily determined by the radial motion of optically 
thick H2O layer. 
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